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The dimension two gluon condensate has been used previously within a simple phenomeno- 
logical model to describe power corrections from available lattice data for the renormal- 
ized Polyakov loop and the heavy quark-antiquark free energy in the deconfined phase of 
QCD P, 01 • The QCD trace anomaly of gluodynamics also shows unequivocal inverse tem- 
perature power corrections which may be encoded as dimension two gluon condensate. We 
analyze lattice data of the trace anomaly and compare with other determinations of the 
condensate from previous references, yielding roughly similar numerical values. 

Introduction. For zero and for infinite quark masses (gluodynamics) QCD is invariant under 
scale and conformal transformations at the classical level. This classical invariance is broken, 
however, by quantum corrections due to the necessary regularization of ultraviolet divergences 
which introduces a mass scale, Aqcd; the divergence of the dilatation current equals the trace of 
the improved energy- momentum tensor 0] yielding the so-called "trace anomaly" At finite 
temperature, the energy density e and the pressure p enter as [1, [1, 0, 

e-3p=^((G;;j2)^(e;i), (i) 

where G^i, ~ d^Ai, — d,yA^ + iglA^, A^] is the field strength tensor and /3(g) — ^dg/dfi is the beta 
function. Far from the conformal limit, where e = 3p, A = (e — 3p)/T'^ is a dimensionless quantity 
providing a measure of the interaction, so it is commonly known as "interaction measure ". A good 
knowledge of A is crucial to understand the deconfinement process, where the non perturbative 
(NP) nature of low energy QCD seems to play a prominent role. In this contribution we analyze 
the highly NP behaviour of the trace anomaly just above the phase transition and describe it in a 
way that is consistent with other thermal observables (see Q for further details). 

Thermal power corrections in gluodynamics. The interaction measure was computed one 
decade ago on the lattice by the Bielefeld group for gluodynamics [13] ■ Fig.[T]shows the lattice data 
for A = (e — 3p)/T^ as a function of T/Tc- A is very small below Tc, because the lightest glueball 
is much heavier than Tc « 270 MeV. It increases suddenly near and above Tc by latent heat of 
deconfinement, and raises a maximum at T « 1.1 Then it has a gradual decrease reaching zero 
in the high temperature limit. The high value of A for Tc < T < (2.5 — 3)Tc corresponds to a 
strongly interacting Quark-Gluon Plasma picture. 

From our previous experience [l|, Q and following a remark by Pisarski ^lli] , in Fig. [T] we plot 
(e — 3p)/T'* as a function of l/T^ (in units of Tc) exhibiting an unmistakable straight line behaviour 
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FIG. 2: The temperature dependence of the renormahzed Polyakov loop in units of the critical temperature. 
Lattice data are from [ij, for A'^^ x Ny — 32^ x 4 and 32^ x 8. We plot the perturbative result at LO and 
NLO, and the fit using apoi + bpoi {Tc/Tf. 

in the region slightly above the critical temperature, of the form 

A = (e - W)/T^ = fltra + fetra [Tc/T f , (2) 

and corresponding to a "power correction" in temperature. A fit of the lattice data [N^ x Nr = 
323 y. 8) foj, y/y^ > 112, yields atra = -0.02(4), 6tra = 3.46(13), x^DOF = 0.35. Power 
corrections also appear in e and P, just by applying the standard thermodynamic relations. 

This behaviour clearly contradicts perturbation theory (PT) which contains no powers but only 
logarithms in the temperature, a feature shared by hard thermal loops and other resummation tech- 
niques (see e.g. [H,[l3]), explaining why they have failed to describe lattice data of the free energy 
below iTc- Fig. [5] shows the lattice data of Polyakov loop from Ref. suggesting again P, |3| a 
linear fit of the form —2 logL = a^ai + ^poi {Tc/Tf . In what follows, we show a phenomenological 
model that describes consistently all these power corrections in an unified way. 

Dimension two gluon condensate. The gluon condensate (G^) = g'^{{G'^^f) describes 
the anomalous (and not spontaneous) breaking of scale invariance, and hence is not an order 



parameter of the phase transition. ActuaUy, the order parameter is the vacuum expectation value 
of the Polyakov loop which signals the breaking of the Z{Nc) discrete symmetry of gluodynamics as 
well as the deconfinement transition. A dimension two gluon condensate naturally appears from a 
computation of the Polyakov loop in which a Gaussian distribution of eigenvalues is considered. In 
the static gauge, 9o^o(x, xq) = 0, this Gaussian-hke, large Nc motivated, approximation gives [l| 



L(r) = /^tr,e'«^«W/^Wexp 



ANrT^ 



(3) 



valid up to 0{g^) in FT. Aq is the gluon field in the (Euclidean) time direction. From here it is 
immediate to relate the Polyakov loop to the gluon propagator in the dimensionally reduced theory 



SabTD, 



3(k) = J d^x{AoA^)AoAy))e 



-ik-(x-y) 



(4) 



The dimension two gluon condensate (7^(Aq ^) is obtained from Eq. ^ in the limit x ^ y. The 



perturbative propagator L)QQ(k) 



being mo ^ T the Debye mass, leads to 



the known perturbative result of Gava and Jengo |15j , which fails to reproduce lattice data below 
6Tc. A NF model is proposed in Ref. [l| to describe the lattice data of the Folyakov loop, and it 
consists in a new piece in the gluon propagator driven by a positive mass dimension parameter: 



(5) 



This ansatz parallels a zero temperature one [16j, where the dimension two condensate provides 
the short-distance NF physics of QCD and at zero temperature this contribution yields the well 
known NF linear term in the qq potential. A justification of Eq. ([5]) based on Schwinger-Dyson 
methods has been given • The new propagator generates a NF contribution to the condensate, 
{Ala) = (Ala)'' + (Ala)'"'', which is related to through (A^^J^p ^ (^2 _ (g^rmz,), 
where rho = mo/T, so that it leads to the thermal power behaviour that we observe in Fig. [21 
The Gaussian approximation has also been used in Ref. to compute the singlet free energy of a 
heavy qq pair [3, through the correlation function of Folyakov loops. 

Non perturbative contribution to the Trace Anomaly. The model of Eq. ([5]) can be 

easily used to compute the trace anomaly Eq. ([T]) in gluodynamics. Assuming the leading NF 
contribution to be encoded in the A^ a field and taking Ai a — yields 



\NP 



(6) 



The r.h.s. is obtained from Eq. ^ by expanding in the limit x ^ y and looking at the quadratic 
term in r = |x — y|. Note that the NF model is formulated in the dimensionally reduced theory, 
so the gluon fields are static. This formula produces the thermal power behaviour of Eq. ^ with 



6traT^ = -3™^(A^,jNP^(g)/5. 



(7) 



If we consider the perturbative value of the beta function (3{g) ~ + 0{g^), the r.h.s. of Eq. ([7|) 
shows a factor in addition to the dimension two gluon condensate g'^{A^ a)^^ ■ the fit of the 
trace anomaly data is sensitive to the value of the smooth T-dependent 5, without jeopardizing the 
power correction. For the Folyakov loop the sensitivity in g is only through the perturbative terms, 
which are much smaller than the NF ones. When we consider the perturbative value gp up to 2- 
loops, we get from the fit of the trace anomaly g'^{A^ a)^^ = (2.63 ± 0.05 Tc)^, which is a factor 1.5 
smaller than what is obtained from other observables. This disagreement could be partly explained 
on the basis of certain ambiguity of 5 in the NF regime. A better fit of the Folyakov loop and heavy 
quark free energy lattice data in the regime Tc < T < ATc is obtained for a slightly smaller g than 



Observable 


g (Ao,a) 


Polyakov loop [1] 
Heavy qq free energy [2] 
Trace Anomaly 


(3.22 ± 0.07 Tc)"" 
(3.33 ± 0.19 Tcf 
(2.86 ± 0.24 Tcf 



TABLE I: Values of the dimension two gluon condensate from a fit of several observables in the deconfined 
phase of gluodynamics: Polyakov loop, singlet free energy of heavy quark-antiquark and trace anomaly. 
Values are in units of Tc- We show the fit for lattice data with Nt = 8. Error in last line takes into account 
an indeterminate value of the coupling constant g — 1.26 — 1.46, being the highest value the perturbative 
gp up to 2- loops at T = 2Tc. The critical temperature in gluodynamics is Tc = 270 ± 2MeV |14i] . 

gp, i.e. g = 1.26 - 1.46 [2]. Taking this value we get from Eq. g^{Al a)^'^ = (2.86 ± 0.24 
in better agreement with other determinations, see Table HI Nonetheless, an alternative method to 
compute the trace anomaly based on the direct computation of e — 3p from the partition function 
of gluodynamics does reproduce the power correction, however with different coefficients f9|. 

Summary and conclusions. The trace anomaly in gluodynamics shows, near and above the 
critical temperature, a clear pattern of power corrections which cannot be matched to perturba- 
tion theory or resummations thereof. It can instead be explained in terms of a dimension two 
gluon condensate whose numerical value agrees with other determinations based on other thermal 
observables and it is also remarkably close to existing studies at T = (see e.g. Refs. [Tol [20|). 
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